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Abstract 

The dynamics of the low energy excitations in a ferromagnet is studied in case a temperature 
gradient is coupled to the local magnetization. Due to the different time scales of changing temper- 
ature and magnetization it is argued that only the coupling between the spatially varying part of 
the temperature field and the magnetization is relevant. Using variational principles the evolution 
equation for the magnetic system is found which is strongly influenced by the local temperature 
profile. The system offers damped spin wave excitations where the strength of damping is de- 
termined by the magneto-thermal coupling. Applying the model to nanowires it is demonstrated 
that the energy spectrum is significantly affected by the boundary conditions as well as the initial 
temperature distribution. In particular, the coupling between temperature and magnetization is 
expected to be several orders stronger for the open as for the isolated wire. 
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I. INTRODUCTION 



Spin-based electronics is a topical and challenging subject in solid-state physics includ- 
ing the manipulation of spin degrees of freedom, in particular the spin dynamics and the 
spin-polarized transport [Ij. In this regard it is desirable to find proper means for the 
creation and control of spin currents. As it was firstly observed in [2j a spin current can 
be thermally generated by placing a magnetic metal in a temperature gradient. This phe- 
nomenon is called spin Seebeck effect (SSE) and can be detected by attaching a Pt wire on 
top of the magnetic metal and measuring a transverse electrical voltage in the Pt wire via 
the inverse Spin Hall effect, the conversion of a spin current into an electric current [31 II]. 
Now a new subfield has emerged called 'spin caloritronics' [5] including the SSE as a local 
thermal transport effect [6j. In particular, the SSE combines spin degrees of freedom with 
caloric properties. Further, the SSE was also observed in ferromagnetic semiconductors [7] 
and insulators [8]. Despite the absence of conduction electrons, a magnetic insulator can 
convert heat flow into a spin voltage. The occurrence of a temperature-gradient induced 
spin current along the magnetization direction in a ferromagnet is called longitudinal SSE 
[S]. Here, the magnon-induced spin current is injected parallel to the temperature gradient 
from a ferromagnet into an attached paramagnetic metal. Theoretically the interplay be- 
tween temperature gradient and spin current within the SSE is investigated intensively. So, 
the SSE is attributed to a temperature difference between the magnon temperature in the 
ferromagnet and the electron temperature in the attached Pt contact [10]. Whereas in the 
conventional charge assisted Seebeck effect the Seebeck coefficient can be expressed by the 
electronic conductivity, a microscopic model for spin transport driven by the temperature 
gradient was proposed in [TT]. The authors find thermally driven contributions to the spin 
current in ferromagnetic metals. In case of ferromagnetic insulators a linear-response theory 
was formulated [12] where the collective spin- wave excitations are taken into account. A 
detailed numerical study of the SSE based on a modified Landau-Lifshitz-Gilbert equation 
is presented in In addition, the scattering of a spin current generated by a temperature 
gradient is analyzed in [Ti] where the setup initiates a spin current torque acting on the 
magnetic scatterer. As pointed out in [15] the domain wall motion can be also originated 
in magnonic spin currents. Likewise it is stressed in [16] a spin-transfer torque in a multi- 
layer structure can be initiated by thermal transport from magnons. Recently [17] reported 
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that the SSE is considerably enhanced by nonequihbrium phonons that drag the low-lying 
spin excitations. Due to a phonon-magnon drag mechanism one finds a phonon driven spin 
distribution [TH] which infiuences the SSE. The amplification of magnons by thermal spin- 
transfer torque is observed in [19] where the system is subjected to a transverse temperature 
gradient. Another aspect of such collective magnetic excitations in a magnetic insulator 
is the transmission of electrical signals |20]. Measurements of the Seebeck coefficient and 
the electrical resistivity of polycrystalline magnetic films are presented in [21]. The results 
provide important information for further studies of the SSE in nanomagnets. Whereas 
most studies of spin caloric effects are obtained for thin films, patterned ferromagnetic films 
demonstrate the profound effect of substrate on the spin-assisted thermal transport [22j. 

Motivated by recent experimental and theoretical findings related to SSE we study the 
magnetization dynamics under the infiuence of a temperature gradient. Because the relevant 
time scale of the magnetization seems to be totally different from the time scale of the tem- 
perature field our analysis is focused on the case that the temperature gradient is preferably 
coupled to spatial variation of the magnetization. Due to spin-wave excitations those spatial 
variations of the local magnetization have an impact on the spin dynamics. To that aim 
we propose an action functional including both fields, the local magnetization field m(x, t) 
as well as the local temperature field T(x, t). The functional includes a symmetry allowed 
coupling between the fields. Using variational principles we derive the relevant Bloch-type 
equations of motion. The evolution equations are discussed in detail for one-dimensional 
structures under distinct boundary conditions. In detail we study a ferromagnetic nanowire. 



II. DERIVATION OF THE EQUATION OF MOTION 

Let us consider a ferromagnet below its Curie temperature under the infiuence of a 
temperature field varying in space and time. Then the action for the total system is composed 
of three parts 

S = Sm + ST + Sl. (1) 

The magnetic part Sm giving rise to the local precession of the magnetic moments m(x, t) 
reads [231 
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Here, the Jap characterize the exchange interaction between the magnetic moments while 
the vector field A is a functional of the magnetic moments m as indicated in Eq. ([2]). It 
will be specified below. Because our system is confined in a finite volume it is appropriate 
to separate the scalar temperature field T(x, t) = $(x)0(t). Following [21] the temperature 
part is written as 

The constant k is the heat conductivity or temperature diffusivity while k is the separation 
parameter with the dimension of a wave vector. Its physical meaning will be discussed 
below. The separation of the temperature field is likewise suggested by the well separated 
time scales of the magnetization and the heat conduction. For that reason we argue that the 
interaction part Sj includes only the coupling of the magnetization to the spatial varying 
part of the temperature field: 

With other words the coupling refiects the assumption that the spatial variation of the 
magnetization is directly coupled to the gradient of the x-depending part $(x). The tensor 
Fo^^ refiects the coupling strength. As noticed above the separation parameter k in Eq. ^ is 
proportional to an inverse length scale k ^ describing the range where the temperature 
is significantly changed in space. Having a nanowire or thin films in mind we assume 
/im < It ^ iTim defined by the length extensions of the sample. The characteristic time 
Tj, for pure heat conductance in a confined medium is ~ {K,k^)~^ ~ Taking the 

thermal diffusivity k = 2.51mm^s~^ for Y3Fe50i2 (YIG crystal) \2E\ one estimates the time 
scale Tt in the range fis < Tj, < s. In contrast magnetodynamics occurs a considerably 
shorter time scale within the nanosecond range: < ns, i.e. is supposed to differ from 
Tm about a few orders of magnitude, ^ Tt- From Eqs. ([l])-(|4]) one finds the equations of 
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(3) 



(4) 



motion by variational principles: 



+ = , (5a) 



i§ + *^e = r„,,^, (5b) 

K at 0Xj30X^ 



dnia dmp J dt '^^'^ dxpd^ ~^ '^ap^y ^^^^^^ • (^^) 
This set of equations describes the influence of a temperature profile on the magnetization, 



Eq. (5c), and the feedback of the magnetization on the temperature, Eq. (5b). In a first 



approximation we neglect the latter one. Consequently, Eqs. (5a) and (5b) represent the 



conventional heat equation after the separation of the variables. In the following we consider 



an isotropic magnet characterized by Jap = J Sap, see Eq. (5c). Multiplying Eq. (5c) with 
m vectorial it results 

m„ (m ■ curljnA) = J (m x V^m)^ + A^^j^^m^ . (6) 

Here there appears a tensor Kap-ys-, which includes the coupling between magnetization and 
temperature profile according to Eq. (|4]). Since A is symmetric with respect to the indices 
7 and 5 one gets for isotropic systems Kaji-^s = Ci6apS^s + C2{5a'y5ps + dasSp-y). If A = 
Eq. ^ should conserve the spin length m2 = 1 because of the absence of dissipative forces. 
From here we conclude that a pure reversible equation of motion can be obtained by setting 
curlmA = —7^^ m where 7 is the gyromagnetic ratio. As the result one gets the Landau- 
Lifshitz equation without relaxation. Using the same ansatz also for nonzero A and setting 
in leading order Ci = —C and C2 = one finds from Eq. (|6|) the following equation 



^""^^'^^ = -7 J m(x, t) X V2m(x, t) + 7 C [V'$(x)] m(x, t) . (7) 

The sign of the coupling parameter C between temperature and magnetization will be dis- 
cussed in detail below. The final evolution equation Eq. ([T]) does not conserve the spin 
length anymore and has the form of a Bloch-Bloembergen-like equation |26l[27j. Therefore 
it should be applicable for ferromagnetic systems on a mesoscopic scale. To affirm the rea- 
sonability of the result in Eq. ([T]) we briefly outline an alternative way of its derivation. Now 
the starting point is the energy functional of the total system 

J- = + J-T + J-/, (8) 
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where it is likewise decomposed into a magnetic part Tm-, a temperature part Tt and an 
interaction part J-}. Within the Ginzburg-Landau theory the functional is given by 

J V 2 dxa dx^ 2 dxadxp "^'^ dxp dxjj 
where the coupling tensors Jafs and Tap-y are the same as before in Eqs. ^ and Q. Moreover, 

we have chosen a quadratic coupling of the gradient of the temperature characterized by 

Xai3- The function $ = $(x) is again the space-dependent part of the scalar temperature 

field T(x, t) = $(x)0(t). In order to derive the desired equations of motion we follow the 

common approach 



Pm J d y {mc,{x),mp{y)} ^^^^y-^ = pmea/s-ym/sj— . (10b) 

Here, the Poisson brackets {mQ,(x), m^(?/)} = 6 {x — y) ea/3^m^{x) for classical magneti- 
zation densities are introduced, see Now we use the separation ansatz for the 
temperature in Eq. (10a) and restrict the calculation to isotropic situations; Jai3 = JSai3, 



= —C 6ai3S-ys and additionally Xai3 = X^aj3- Identifying further xPt = and pm = —7 
the resulting evolution equations coincide exactly with Eqs. (5a), (5b) and ^ for the tem- 



perature and the magnetization, respectively, as derived from the expressions for the action 
in Eqs. ([T])-(|4]). The advantage of the second method is that the firstly unknown quan- 
tity A(m), introduced in Eq. ([2]), is not necessary. On the other hand, the coincidence of 
the obtained results for the equations of motion based on different approaches justifies the 
assumptions made for the vector field A. 

III. APPLICATION TO FERROMAGNETIC NANOWIRES 

In this section Eq. ([T]) is analyzed for a ferromagnetic nanowire as schematically de- 
picted in Fig. [1} Because the model is focused on the ferromagnetic phase and the related 
collective low energetic excitations, the magnetization-field is decomposed into m{x,t) = 
mo + il}{x,t) = ipi{x,t)ei + ip2{x,t)e2 + mo 63. Here ttiq is the constant magnetization 
whereas the spin- wave amplitudes ipii'4'2 obey the coupled equations 

Ot OX^ OX'^ 

= J mo —— + 7 C — - iIj2 . 
at ox^ OX'' 
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FIG. 1. (Color online) Nanowire with locally varying temperature. 
The magnetization dynamics is accordingly influenced only by the spatially varying temper- 



ature fleld $(a;), the complete dynamics of which are given by Eqs. (5a) and Eq. (5b). In 
one dimension it results 



ox'^ at 



(12) 



The solution for $(a;) in Eq. (12) can be written as a superposition $(x) = ^„ ^n{x). Hence 



d'^^/dx^ = — J2n ^n^n cau bc substituted in Eq. (11). In terms of the Fourier transformed 



^n{(l^t) = •^7~{$(x,t)} and (pi{q,t) = J-'T{'4>i{x,t)} Eq. (11) can then be transformed into 
d 



dt 

d_ 
dt 



(pi{q,t) = -a q (p2{q,t)-ra 



ip2{q,t) = a q (pi{q,t) - r a' 



/ dq'Mq-q',t)Mq: 
X] j dq' y^2{q - q',t) <l>„(g'; 



(13) 



Here the time t is measured in terms of t = 7 J mot, where for convenience t is already 



called t in Eq. (13). The other dimensionless quantities are deflned by 

J = J , C = K niQ , r = K/ J , 



(14) 



where a is the lattice constant. The separation parameters kn and the ^n{q) are determined 



as the solution of Eq. (12) under consideration of the relevant boundary conditions. 
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A. Zero heat flow at the boundaries 



Let us firstly study the case when there is no heat current flowing through the edges of 



the wires. The initial and boundary conditions for the heat equation Eq. (12) read 

dT{x,t) 



T(x,t = O) = 0o5^$. 



dx 



0. 



x=0,L 



(15) 



where L is the length of the nanowire. The solution of Eq. (12) fulfllling the boundary 



condition Eq. (15) is 

<l>„(x) = Un COs{knX) , kn = 

Performing the Fourier transformation we obtain 



n vr 



n = 0,1,2,.... 



$n(g) = Mn TT [Hi + ^n) + Hi - ^n)] ■ 



(16) 



(17) 



Inserting this relation into Eq. (13) leads to terms of the form (pi(q ± kn) which can be 



expanded with respect to the small parameters kn- 

dipi 



ip,{q ± kn) ~ V>iil) ±-^kn + {kl) . 



(18) 



Is A the wave length of the low energetic excitation in the ferromagnet with the corresponding 



wave vector g ~ 1/A and setting fc„ ~ 1/L the expansion made in Eq. (18) is valid for 



/cn^g^l/aorL^A^a, respectively. Introducing t] = ipi + iip2 and combining 



Eqs. (17)-(18) with Eq. (13) one obtains the following evolution equation for rj: 



d 



dt 



vVil, ^ = (i a^?^ - ") Vil, t), a = 2nra'^ ^ k^ 



Ur 



The solution for the Lp = {ipi,ip2) with the initial conditions ip{t = 0) = (p^ = (v5oi)V^02) is 
written in terms of the real time t as 

I cos [u (q) t] -sin [uiq)t]\ ,22 n 

(p{q,t)= cfiQexp[-at\, u{q) = J moa q . (19) 

\ sin [uj{q) t] cos [uj{q) t] I 

Here the spin wave energy uj{q) and the damping constant a are deflned by 

uj{q) = 1 J moq"^ , a = 27frJ'jmo^klun ■ (20) 

n 

The expression for uj{q) reflects the common quadratic spin wave dispersion relation. The 
temperature gradient leads to damped spin waves whereas the damping parameter a is 
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FIG. 2. (Color online) Spin waves (fi{q,t)/(fo (color scale) for a = lA, L = 1/xm and r = 100 
for the cases: (a) vanishing heat current at the edges of the nanowire according to the initial heat 



profile I in the inset in Fig. 3(a) , (b) fixed temperatures at the edges referring to the initial heat 



profile III in the inset in Fig. 3(b) 



defined in Eq. (20). The strength of the attenuation a can be estimated from Eq. (20) 
and Eq. (14) using a oc C*^^fc^M„. In order that the damping is physically realized we 



have to discuss the sign of the parameter C oc r introduced in (j?]). It characterizes the 
coupling between the local magnetization and the local temperature. Moreover, the sign of 
Yin '^n IS determined by the sign of the coefficients m„ which is given in terms of the initial 



temperature distribution, compare Eqs. (15)-(16). A positive damping parameter a > is 
realized by assigning 



r ^ 0. 



(21) 



The t- and g-dependent ratio (pi/(po is plotted in Fig. 2(a) As expected the frequency is 
increased when q is enlarged. Moreover the damping should be observed within a picoseconds 
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FIG. 3. (Color online) Spin wave solutions for a = lA, L = l^m and g = 4 x lO'^m^^ for the cases: 
(a) vanishing heat current at the edges of the nanowire with the initial heat profiles I (n„, = 
uq = 10) and II (un = U4 = —1, uo = 10) mentioned in the text; (b) fixed temperatures at the 
edges with the initial heat profiles III {vn = — 1/n^, Ti/Gq = 12, T2/0O = 8.5) and IV {vn = vs, 
Ti/Qq = 10 = T2/Q0) in the inset mentioned in the text. Here only the envelope of the solution 
referring to heat profile III is plotted. 



time-scale. To illustrate that we present the dynamic solutions for different heat profiles in 



Fig. 3(a), The initial temperature distribution is chosen by Tq/Qq = uq + Xln""" cos[/c„x] 



and is depicted in the inset of Fig. 3(a) To reach a damping of the spin wave in the ps- 
range we set |r| = K/J = 100, where r is a measure for the relation of the coupling of the 
temperature and magnetization profiles (i.e. V$Vm) and the magnetization-magnetization 



(i.e. VmVm) coupling, see Eq. (14). In that case one observes a strong coupling between 
temperature and magnetization. 
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B. Fixed temperatures at the boundaries 



As the second realization we consider the case that the edges of the wire are held at fixed 
temperatures: 



T{x,t = 0)= eo^$„(x), $(x = 0) = $(x = L) = 0. 



(22) 



Therefore, Eq. (22) refers to the transformed temperature field T{x,t) = T{x,t) — (1 — 



x/L) Ti — X/LT2, where Ti = T{x = 0, t) and T2 = T{x = L, t) are fixed and T(x, t) is the 
original temperature profile in the wire. The spatial temperature profile is obtained 



from Eq. (12) combined with the boundary conditions in Eq. (22). The solution is the 



superposition of 



n IT 



^n{x) = Vn Sm{knX) , ^ ' ^ = ^'2, .... 



(23) 



The Fourier transform of that $„(a;) reads 



= it;„vr [6 {q + kn) - S {q - A;„)] . 



(24) 



Utilizing again the expansion in Eq. (18) and the transformation r] = ipi + i{p2 v/e find by 



inserting Eq. (24) in Eq. (13) 



d - d 



(25) 



This partial differential equation of first order can be solved by using the method of char- 
acteristics. The solution for if = {(pi,(f2) with the initial conditions ifi{t = 0) = c^q is given 
by 



cos[il{q,t)t] — sm[il{q,t)t] 
(p{q,t) = I I (Pq exp[J'jmoP{t)qt]. 

sin [r2(g, t) t] cos[fi(g,t)t] 



(26) 



The quadratic spin wave dispersion relation fi(g, t) and the damping parameter (3{t) are 
defined by 



fi(g,t) = J'jmo ( q 



Pit) = 27rrJ7mot ^klvn- 



(27) 



Obviously the quadratic spin wave dispersion relation is now modified by a time depending 
correction term JjmoP'^{t)/3, i.e. the spin wave dispersion is changed by the damping 
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parameter P{t) which occurs as a factor in the exponential function in Eq. (26). The decay 
of this function is guaranteed for negative values /3(t) < 0. In the same manner as before, 



using Eq. (27) and Eq. (14), we estimate /3 oc C'Xln^n'^n- Thus, one has to distinguish 



r ^ 0. 



(28) 



Referring to the initial heat distribution (T was introduced below in Eq. (22)) we assume 
Tq/Qq = (1 ~ x/ L)Ti/Qq + XT2/L60 + J2'^n sin[A;„x], see the inset in Fig. 3(b) The 



solution (pi(q,t)/(po is shown in Fig. 2(b), Contrary to the previous case of an absent heat 



flow through the ends of the wire now we observe an enhancement of the lifetime of the spin 
waves, cf. the time scales of Figs. |2(a)| and 2(b)[ The same effect of increased spin wave 
lifetimes is noticeable in Fig.[3]where the solution for q = 4x 10®m~^ is illustrated for different 
heat profiles but identical values of the coupling relation r for both types of boundary 



conditions. In Fig. 3(b) we only plot the envelope of the solution for one of both cases. 



However, the curve shown is featured by a decrease of the frequency up to a time of ~ 430ps 
which is determined by the maximum of the argument Q{q,t)t of the periodic functions. 
Thereafter the frequency increases continuously with time. A quantitative comparison of 
the lifetimes for comparable temperature distributions in Fig. |3] offers that the characteristic 
time scale is 40 — 90 times larger for the case when heat flows through the edges than for 
the case the nanowire is isolated. The reason for that consists of the expansion made in 



Eq. (18). Because in Eq. (25) all even terms cancel mutually in the combination (pi(q + 



kn) — <^i{q — kn) and hence only the odd terms remain. This leads to an alteration of 



the damping constant in Eqs. (26)-(27) compared to a in Eqs. (19)-(20). The physical 



picture behind might be as follows: As discussed the type of the temperature profile $(x) 
determines the strength of coupling. In our model the solution for $(x) is selected by the 



boundary conditions either in Eq. (15) or in Eq. (22). The complete isolation of the nanowire 



imposes the restriction on the thermo-magnetic system that it is not able to exchange energy 
with the environment. This leads to very limited possibilities for the self-organization of the 
magneto-caloric behavior. Otherwise, if heat is allowed to flow through the edges of the wire 
the thermo-magnetic system is in contact with its surrounding. Thus the self-organization 
process can be optimized due to a broader variety of its realization. 
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IV. CONCLUSIONS 



In the present work the low energetic excitations in a ferromagnetic nanowire had been 
analyzed when the spatial variations of the magnetization is coupled to the correspond- 
ing temperature gradient. The resulting spin- waves are damped due to the temperature 
gradient. This result contributes to the recent effort in understanding spin caloritronics 
phenomena. Whereas the spin Seebeck effect is focused on the creation of a spin current 
due to a temperature gradient the goal of our study is to elucidate the behavior of the per- 
manent effect of that temperature gradient on the excitation of the magnetic system. Once 
the magnetic system is excited the shape of the temperature profile influences the damping, 
too. Further, the external conditions under which an experiment is realized should affect 
the dynamics as well. In this regard we found that the system is strongly controlled by the 
boundary conditions which modify the damping parameter. While the complete isolated 
nanowire offers conventional spin wave excitations which are damped due to the coupling to 
a spatial changing temperature field the situation is more complex in case the ends of the 
nanowire are held at fixed temperatures. Referring to the latter realization it results that 
the damping parameter is time-dependent and modifies the spin wave dispersion relation. 
However, if one considers systems which are subjected to identical boundary conditions the 
actual initial temperature distribution also affects the spin dynamics. In principle it should 
be possible to validate our theoretical findings experimentally by spin Seebeck effect mea- 
surements or ferromagnetic resonance techniques in low dimensional samples. Our studies 
can be accordingly extended to higher dimensions. Further, the model can be refined by 



taking into account higher order terms in the expansion made in Eq. (18). The inclusion of 
stochastic forces should also improve the model slightly. 

One of us (T.B.) is grateful to the Research Network 'Nanostructured Materials' , which 
is supported by the Saxony-Anhalt State, Germany. 
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